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Abstract 

By a discrete torus we mean the Cayley graph associated to a finite product of finite cycle 
groups with generating set given by choosing a generator for each cyclic factor. In this article 
we study the spectral theory of the combinatorial Laplacian for sequences of discrete tori 
when the orders of the cyclic factors tend to infinity at comparable rates. First we show 
that the sequence of heat kernels corresponding to the degenerating family converges, after 
re-scaling, to the heat kernel on an associated real torus. We then establish an asymptotic 
expansion, in the degeneration parameter, of the determinant of the combinatorial Laplacian. 
The zeta-regularized determinant of the Laplacian of the limiting real torus appears as the 
constant term in this expansion. On the other hand, using a classical theorem by Kirchhoff 
the determinant of the combinatorial Laplacian of a finite graph divided by the number of 
vertices equals the number of spanning trees, called the complexity, of the graph. As a result, 
we establish a precise connection between the complexity of the Cayley graphs of finite abelian 
groups and heights of real tori. It is also known that spectral determinants on discrete tori 
can be expressed using trigonometric functions and that spectral determinants on real tori 
can be expressed using modular forms on general linear groups. Another interpretation of our 
analysis is thus to establish a link between limiting values of certain products of trigonometric 
functions and modular forms. The heat kernel analysis which we employ uses a careful study 
of J-Bessel functions. Our methods extend to prove the asymptotic behavior of other spectral 
invariants through degeneration, such as special values of spectral zeta functions and Epstein- 
Hurwitz type zeta functions. 



1 Introduction 



1.1. The problems we study in the present article begin with the following very elementary 
question. For any d > 1, let N = (ni, • • • , i%d) denote a d-tuple of positive integers, and consider 
the product 

D(N) = Yl {2d-2cos(2nk 1 /n 1 ) 2cos(2nk d /n d )) ; (1) 

where the product is over all d-tuples K = (fc 1; • • • , kd) of non-negative integers with kj < rij, 
omitting the zero vector in the product. The basic question is the following: What is the asymp- 
totic behavior of D{N) as N -> oo? Setting V{N) = n 1 ---n d , then (log D(N))/V(N) can be 
interpreted as an improper Ricmann sum, and we have the limiting formula 

2 cos(27rx,j)) dxi ■ ■ ■ dxd as each rij — > oo. 

(2) 

*We are grateful to the HausdorfT Center in Bonn for the opportunity for us to meet in the summer of 2007. The 
first-named and second-named authors acknowledge support from NSF and PSC-CUNY grants. The last-mentioned 
author thanks VR, the Goran Gustafsson Foundation, and G.S. Magnussons fond for support. 



V(N) 



log D{N) 



log (2d — 2 cos(27ra;i) 



Z d \R d 
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The integral in ^ exists in the sense of improper integrals from elementary calculus, and the 
convergence of D(N)/V(N) to the improper integral in (HJ) can be verified easily the monotonicity 
of \ogx for x > and calculus. 

One can view D(N) as a determinant of a naturally defined matrix from graph theory. Quite 
generally, associated to any finite graph, there is a discrete Laplacian which acts on the finite 
dimensional space of complex valued functions whose domain of definition is the space of vertices 
of the graph. With our normalization of the Laplacian, defined in section HOI below. D(N) is equal 
to the product of the non-zero eigenvalues of the Laplacian associated to a graph which we call a 
discrete torus. 

The study of D(N) takes on an entirely new level of significance beginning with the 1847 paper 
of Kirchhoff [18J which further recognizes D(N) as a fundamental invariant from graph theory. 
A spanning tree of a graph is a sub-graph which contains precisely one path between each pair 
of vertices of the original graph. Kirchhoff 's Theorem states that the number of spanning trees 
is equal to D(N)/V(N), the product of the non-zero eigenvalues of the Laplacian divided by the 
number of vertices. 

Modern mathematics, theoretical computer science, and statistical physics contain numerous stud- 
ies which in some way involves the number of spanning trees of a given graph, determinants of 
Laplacians, or other symmetric functions of eigenvalues which form basic invariants. With all 
this, we see that the invariant D(N) has considerable significance far beyond the elementary 
considerations which allow for its definition. 

1.2. Summary of the main results. Our work begins with the results from [TS] which estab- 
lishes a thcta inversion formula for the discrete Laplacian acting on the space of vertices of nZ\Z. 
Let I x (t) be the classical /-Besscl function, reviewed in detail in section |2~T1 below . Then using the 
general concept of theta functions with inversion formulas constructed from heat kernels (sec, for 
example, [14j). the authors in |15j prove that for any t > and integer x, the following identity 
holds for the theta function associated to the discrete torus nZ\Z: 

1 n — 1 oc 

n (t,x) = i^ e -(a-ac«(a»fc/n))tf!hrifcx/ n ■ = e -» ^ I x+jn (2t), (3) 

k— j — — oo 

where I x (t) denotes the /-Bessel function, see section 12.21 The generalization of ([3]) to the d- 
dimensional discrete torus, defined as the product space 

d 

DT N = Y[ nj Z\Z, 

3=1 

comes from taking a c?-fold product of the theta functions in By computing an integral 

transform of the (i-dimcnsional theta inversion, we obtain the following formula, given in Theorem 
13.61 For any s £ C with Re(s 2 ) > 0, we have that 

J2 log (s 2 + Aj) = V(N)l d (s) + H N (s) (4) 

where 

{Aj} = {2d — 2 cos(27rfc 1 /ri 1 ) — • • • — 2 cos(2tt k^/n^)} with kj = 0, ■ ■■ , rij — 1, 

as in JTJ, 

oo 

Us) = J (e- 2dt e^ 2 %(2^- e -*) * 
o 

and 

oo 

H N (s) = -J (V s2 < [9 N (t) - V(N)er 2dt I Q (2t) d - 1] + e- f ) j. 
o 
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We now consider a sequence of integral vectors N(u) = (rii(u), ...,nd(u)) parameterized by u € Z 
such that rij(u)/u — * otj as u — * oo for each j. We let ^4 be the diagonal matrix with the numbers 
a, on the diagonal and let V(A) = a% ■ ■ ■ a<j ^ 0. Through a careful study of the infinite series of 
7-Bcsscl functions in ([3]), we obtain the following theorem, which is one of the main results of the 
present article. 

Main Theorem. Let logdet* &dt.n(u) be the log- determinant of the Laplacian of non-zero 
eigenvalues on the d- dimensional discrete torus associated to N(u), and logdet* A.rt.a be the 
log- determinant of the Laplacian on the real torus AZ d \R d . Then 

logdet*A jDT:JV( „ ) = V(N(u))X d (0) + logw 2 + logdet* A rt ,a + o(l) asw^oo (5) 

where 

(oo 
\og2d- J (e- 2dt (I (2t) d - 
o 

There are a number of facets of ([3]) which are interesting. First, the determinant det* A^j^jv^) 
is a legitimate finite product of eigenvalues, whereas det* Art, a is defined through zeta function 
regularization. From this point of view, we have connected a zeta regularized determinant with 
a classical determinant. Moreover, since the first terms in ([5|) are universal, this allows for the 
possibility of transfering knowledge, for example, between the minimal spectral determinants for 
real tori, as studied in [51 123j. to minimal spectral determinants for discrete tori, or alternatively, 
the minimal number of spanning trees. In this context it is relevant to remark that the main 
theorem holds without any changes for general tori, see section 17.51 Second, the discrete torus 
DTjy can be viewed as a lattice on the real torus AZ d \R d , and the degeneration which occurs when 
u tends to infinity amounts to considering a family of discrete tori which are becoming uniformly 
dense in AZ d \R d . The main theorem above proves a type of "re-scaled continuity" when studying 
the asymptotic behavior of spectral determinants. Finally, the classical Kronecker limit formula 
for Epstein zeta functions amounts to the evaluation of logdet* Art, a in terms of a generalization 
of Dcdekind's eta function to a GL(d, Z) modular form. With this, we have established a precise 
connection between the asymptotics of the number of spanning trees on families of discrete tori 
and modular forms. 

The verification of the main theorem in the case d = 1 can be carried out directly; the details are 
presented in section (|6.1|) . When d = 2, additional explicit computations are possible, in which 
case the main theorem becomes the following result. 

Main Theorem in the case d = 2. Let N(u) = (rii(u), 112(11)) = (72,1,712), and assume that 
N(u)/u — > (ai,Q!2). With the above notation, we have, for any integer K > 3 ; the asymptotic 
formula 

AG K ~ 3 
logdet*A £ , TiA r (u ) = n x n 2 V log(nin 2 ) + logfl r?(ia 2 / 'a 1 )\ 4 a 2 / 'at) + 2J F k (u) + 0(u K ) (6) 

fc=i 

as u — > 00, where G is the classical Catalan constant, r\ denotes Dedekind's eta function, and the 
functions Fh(u) are explicitly computable and satisfy the asymptotic bound Fk(u) = 0(u~ k ~ 2 ) as 
u — > 00. 

The Main Theorem (J5j comes from studying the special value s = in J4]). The analysis we 
develop in the proof of (J5j extends to prove the asymptotic behavior of (|4]) for all s G C with 
Re(s 2 ) > as well as for special values of the spectral zeta function J^ . Aj w for w € C with 
Re(u>) > 0. In subsection 17.31 we prove the following theorem. 

Theorem. With notation as above, let £4 be the spectral zeta function on the real torus AZ d \R d . 
Then for any w € C with Ke(w) > 0, we have that 



lim u- 2w 



u— »oo 



/ 

( \ 

~ T{w) J v " J " Y\ -^^J (4Tr) d / 2 (w~d/2)T(w) 

V 



C^M-S Ae-»(V2t)))Vg \=CaH- ,^wZ ( _Ltv...v ( 7 ) 
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In particular, for any w £ C with He(w) > d/2, we have 



u 



lim {u- 2w ( N{u) (w))=( A {w). 



(8) 



1.3. Comparisons with known results. The case d = 2 has been studied since at least 
the 1960s by physicists, starting by Kasteleyn's celebrated computation [16] of the lead term 
AG /it. The next order terms in the asymptotic expansion was later obtained by Barber, see 
Duplantier-David for references and [HI (3.18)] for the precise statements. Going further, in 
[9] the authors establish the error term Olu" 1 ), while our explicit error term expansion begins a 
term of order 0(u~ 3 ). The authors in [3] prove their results by studying the asymptotic behavior 
of the spectrum of the Laplacian, as opposed to the asymptotic behavior of the heat kernel and its 
integral transforms, which is the approach taken in the present article. Subsequent authors have 
studied asymptotics of the Laplacian on more general subgraphs of Z 2 , see, for example, Burton- 
Pcmantle [2] and Kenyon [17]. For general d > 2, the articles [2] [4] [TTJ [24] [25] have discussions 
which give the lead term asymptotics in ([5]). However, there is no discussion in these papers or 
elsewhere regarding the next order term in the asymptotic expansion in ([5]) . 

As stated, in sections 17.31 and 17.41 we apply our analysis to study the asymptotic behavior of the 
families of spectral zeta functions and Epstein-Hurwitz zeta functions through degeneration. These 
problems are considered in [9] in the case d = 2 for general arguments for the Epstein-Hurwitz 
zeta function and for the special values w = 1 and w = 2 of the spectral zeta function (see section 
17.41 for a clarification of the notation). Specifically, equation (3.24) of [3] asserts that, "after a 
rather long algebra" (quoting [3]), one has the asymptotic formula 



as u — ► oo. For comparison, our result ([8]) specialized to d = w — 2, yields the asymptotic formula 



as u — ► oo. Going further, the current literature contains only a handful of other considerations 
beyond ©, all in the case when d — 2, whereas our result is completely general. 

As an aside, note that by equating the lead terms in ([9]) and (|10[) and setting y = oixjoii, we arrive 
at the formula 



A direct proof of pip can be obtained by taking the Fourier expansion of the non-holomorphic 
Eisenstein series E(z, s) for SL(2, Z) with z — iy and s = 2, using the evaluation of the A"-Bcsscl 
function K s _\ji in terms of exponential functions when s = 2. We thank Professor Cormac 
O'Sullivan for clarifying this point for us, and we refer the interested reader to his forthcoming 
article ([H]) for a systematic presentation of identities of this form. 

In Riemannian geometry, the determinant of the Laplacian obtained through zeta function reg- 
ularization has been used extensively as a height function on moduli space; see, for example, [5] 
and [20] . We view the asymptotic expansion in ([5]) as establishing a precise connection to the well 
established notion of complexity from graph theory. Compare with the discussion in |22[ p. 619] 
or [17| p. 242]. Another context where the number of spanning trees in discrete tori appears is 






(10) 




(11) 
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in the study of sandpile models |3J. Going beyond the setting of discrete and real tori, we are 
optimistic that the point of view taken in the present paper will extend to address problems in 
the study of other families of Caylcy graphs of discrete, possibly infinite, groups. 

1.4. Outline of the paper. In section[5J we establish notation and present background material 
from elsewhere in the mathematical literature. 

In section [3J we study the theta function ([3J associated to the action of the discrete Laplacian 

d 

on functions on DTn = Yl «jZ\Z. Following the results from [15j . we begin with the inversion 
j'=i 

formula ^ obtained by expressing the heat kernel on Z d in terms of i-Bessel functions and then 
pcriodize to obtain an expression for the heat kernel on the discrete torus. The main result in 
section [3] is Theorem 13.61 stated in (0} above, which computes the Gauss transform (Laplace 
transform with a quadratic change of variables) of the heat kernel on DTjq. By taking a special 
value of the Gauss transform, we obtain an expression for the determinant of the Laplacian on 
DTpf in terms of integral transforms of I-Bessel functions. 

In section [4j we prove general bounds for i-Bcsscl functions, building from the fundamental esti- 
mates proved in [3T]. Many of the computations in section 0] involve re-scaled I-Bessel functions. 
For instance, we establish asymptotic behavior (Proposition I4.7[) and uniform bounds (Lemma 
14. 6p in the parameter u of 

f(u;x,t) =ue- u2t I ux (u 2 t) (12) 

In particular, we prove that as u tends to infinity, then f(u; x, t) approaches the heat kernel on R 
associated to a certain scaling of the usual Laplacian. 

In section [S]we define and study spectral asymptotics on degenerating sequences of discrete tori. 
Theorem l5.81 and its reformulation in Remark l5.91 proves the main result of this article as stated in 
([5]) , which is the asymptotic behavior of the determinant of the Laplacian on a degenerating family 
of discrete tori. The two key steps are to first employ the change of variable t — ► u 2 t and then, 
afterwards, utilize a careful decomposition of the Bessel integrals involved. The various integrals 
are then analyzed individually and matched up with terms in the meromorphic continuation of 
the spectral determinant on the continuous side. 

Whereas we highlight Thcorcm l5.8l as the main result of the present article, there arc many other 
applications of the analysis from sections [3] and 0] going beyond the study of the determinant of 
the Laplacian which is the focus of attention in section [5] 

In section [SJ we deduce our Main Theorem for d = 2 and show that it agrees with the work in 
[S]. Moreover, we show how to obtain an explicit expansion of the error term. This docs not 
appear in previous works and could be of importance when comparing heights of discrete tori and 
corresponing continuous tori. 

The main points wc address in section[7]are the following: Examination of the numerical evaluation 
for the lead term in various asymptotic expansions for arbitrary dimension; determination of 
asymptotic behavior of spectral zeta functions and Epstein-Hurwitz zeta functions for general 
dimension; and investigation of the notational changes needed to consider general sequences of 
degenerating discrete tori. 

2 Preliminary material 

The purpose of this section is to establish notation and recall relevant material from elsewhere in 
the mathematical literature. 

2.1. Basic notation. We use t £ R + to denote a positive real variable and x S Z to denote an 
integer variable. For any integer n > 1, we call the quotient space nZ\Z a discrete circle, and a 
finite product of discrete circles will be called a discrete torus. The product of d discrete circles 
formed with the integers ni,-- - , will be denoted by DTn where N — (n^j—i,...^- 
For any function / : Z — > R, wc define the Laplacian Az by 



A z f(x) = 2f(x) - (f(x + 1) + f(x - 1)) 
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The Laplacian on Z d is the sum of d Laplacians Az, one for each coordinate. The spectrum of the 
Laplacian Az acting on function on nZ\Z is easily computable (see, for example, |15j); hence the 
spectrum on DTn is simply the set of sums of eigenvalues for each discrete circle. Specifically, let 
{Aj} denote the set of eigenvalues of A acting on function on DT^. Then, with our normalization 
of the Laplacian, 

{Aj} = {2d- 2cos(27rfci/ni) 2 cos(27rfcd/n d )} 

with ki = 0, • ■ • , ni — 1, ■ • • , kd = 0, • • ■ , rid — 1. Let V(N) — n\ ■ ■ ■ rid, which can be viewed as 
a volume of DT^.d- There are V(N) eigenvalues of the Laplacian on DTn, and V(N) — 1 of the 
eigenvalues are non-zero. 

The discrete torus DTn gives rise to a graph by inserting an edge between two points which have 
a single component that differs by one (i.e., nearest neighbor). This is the Cayley graph of the 

d 

group Yi rijZ\7i with respect to the generators corresponding to the standard basis vectors of 
j'=i 

Zi d . A spanning tree is a subgraph such that given any two vertices in DTn, there is precisely one 
path within the subgraph that connects the two vertices. The Matrix- Tree theorem (see e.g. [12] ) 
asserts that 

# of spanning trees = JJ Aj. 

We will use the notation 

det*ADT,Jv = J| Aj 

where dct* A.dt,n denotes the determinant of the Laplacian on the discrete torus DTn omitting 
the zero eigenvalue. 

2.2. The i-Bessel function. Classically, the J-Bessel function I x (t) is defined as a solution to 
the differential equation 

9 d 2 w dw , 9 9 . 

For integer values of x, it is immediately shown that I x = I_ x and, for positive integer values of 
x, we have the series representation 

I (t) = V {t/2)2n+X (13) 

71 — ' 

as well as the integral representation 

7T 

I x (t) = I J e tcos{e) cos{ex)d9. (14) 
o 

The mathematical literature contains a vast number of articles and monographs which study the 
many fascinating properties and manifestations of the /-Bcsscl functions, as well as other Bessel 
functions. As demonstrated in the analysis in [15) . basic to our considerations is the relation 

I x+ x{t)+I x -i{t) = 2j t I x (t), 

which easily can be derived from the integral representation and simple trigonometric identities. 

2.3. Universal bounds for the /-Bessel function. The elementary definition of the /-Bessel 
function leads to a number of precise expressions for I x (t) , two of which are stated in section 12.21 
Unfortunately, the explicit identities do not easily lead to viable estimates for I x (t). Beginning 
with the differential equation which characterizes I x {t), the author in |21j was able to derive very 
precise upper and lower bounds for the /-Bessel function, which we now state. Let 



g x (t) = vV +t 2 ) +x\og 



x + x /(x 2 +t 2 ) 
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Then for alH > and x > 2 we have that 

e-VCVC^ 5 )) < 7(2^) • (x 2 + ^jV^^e-s-W < e i/(V^)). (15) 

These bounds will play an important role in this article. Indeed, in section 4, we will study the 
function g x (t) — t in order to obtain elementary bounds for ^/te - * I x (t) which will be vital when 
establishing our main result. 

2.4. Heat kernels and zeta functions. In [13], the authors established a theta inversion formula 
obtained from studying the heat kernel associated to the Laplacian on nZ\Z. Let Kz(t,x) be a 
function on Z such that 

d 

A z K z {t,x) + —Kz(t,x) = 

and 

lira K%(t, x) = 1 if x — 0, and equal to zero if x 7^ 0. 

General results from [7] and \§\ prove the existence and uniqueness of Kz(t, x) and, as shown, for 
example in |15j . we have that 

K z (t,x) = e~ 2t T x (2t). 

The existence and uniqueness theorems for heat kernels on graphs apply to heat kernels on nZ\Z. 
Bounds for the 7-Bessel function from section 12.31 imply that the series 



J2 e- 2t I kn+x (2t) (16) 



k— — oo 

converges, and hence is equal to the heat kernel on nZ\Z, which we denote by K n z\z(t, x). Denote 
the eigenvalues on nZ\Z and corresponding eigcnfunctions by {A nj } and {4> n ,j} respectively. This 
leads to the identity 

n— 1 00 

^e- tA "^„ J (x)0„ J (O) = £ e- 2t I kn+x {2t). (17) 

3=0 fe=— 00 

In the case n = 1, this yields a classical identity which expresses an infinite sum of 7-Bessel 
functions as an exponential function. The heat kernel on DTn and corresponding theta inversion 
formula is obtained by taking the d-fold product of (fl~7|) . 

2.5. Spectral analysis on real tori. Given a positive definite d x d matrix A, let RTa denote 
the real torus j4Z d \R d . On R we use the variable y to denote the standard global coordinate, 
from which we have the Laplacian 

A R f(y) = ~m 

and corresponding heat kernel Kn(t,y) which is uniquely characterized by the conditions 

A R K R (t,y) + ^- t K R {t,y) =0 



and 



J K n (t, x - y)f(y)dy -> f(x) as f -> 

R 



for any smooth, real- valued function / on R with compact support. The set of eigenvalues of the 
Laplacian on RTa is given by {(2ir) 2 t mA*m} for m € Z d , where A* is the dual lattice to A. The 
eigcnfunctions of the Laplacian are expressible as exponential functions. The theta function 

Q A (t) = e-^- tmA ' m -\ 
mez d 



8 



with t > 0, is the trace of the heat kernel on RT A . The asymptotic behavior of the theta function 
is well-known, namely that 

Q A (t) = V{A)(ATtt)- d /' 2 + 0(e~ c/t ) for some c> as t -► 0, 

and 

Q A (t) = 1 + 0(e~ ct ) for some c > as t — > oo. 

2.6. Regularized determinants on real tori. In very general circumstances, the spectral zeta 
function is defined as the Mellin transform M of the theta function formed with the non-zero 
eigenvalues. Specifically, for seC with Re(s) > d/2, we define the spectral zeta function 

oo 

U(s) = MO A (s) = j^y / (6a (t) - 1) t s f ■ (18) 

o 

The integral in (JT5J) converges for Re(s) > d/2. Let us write 

l 

,dt 



Us) = j^y / (e^(t) - y(A)(4^)- d / 2 ) t s y 

o 
i 

"(V(A)(47rf)- d / 2 -l)t s * (19) 



r( s ) 







r (Q A (t)-l)t° di 



T(s) J " ' t 

l 

Going further, we can carry out the integral in the second term to get 

l 



Us) = j^y / (e^(t) - y(A)(4^)- d / 2 ) i s f 
o 



(s-d/2)T(s) r(s + l) 

^/(e,( t )-i)4 



r( 

1 

By virtue of the asymptotic behavior of Q A (t) as t — > and i — > oo, the expression (|20|) provides 
a meromorphic continuation of Ot(s) to s 6 C. Therefore, using (]20[) . we can study the behavior 
of Ca(s) near s = 0. The integrals are holomorphic near s = 0, and 1/T(s) = s + 0(s 2 ) near s = 0, 
so then we have that Ca(0) = — 1, which is a point that will be used later. In particular, we have 
that 



Ca(0) = / (6x(t) - V(A)(47rt)- d /i) j + T'(l) 







dt 

T 

(21) 

-l V {A){A^ + J (ex(t)-l)* 
1 

which will play an important role in the proof of our Main Result. 

2.7. Kronecker's limit formula. The special value Ca(^) nas a unique place in classical analytic 
number theory, and the evaluation of Ca(0) m terms of modular forms is generally referred to as 
Kronecker's limit formula, which we now describe. The content of this section has its origins in 
[10], and the discussion we present here comes directly from [8]. 
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Let Q denote a d x d positive definite matrix, and u G R d any vector. If we write Q = {qij), then 
the quadratic form associated to Q is denned by 

Q(u) = '/■.."/";• 



The Epstein zeta function associated to Q is defined for Rc(s) > d/2 by the convergent series 

z( S ,q)= Q( m r s - 

mez d \{o} 

From the discussion in sections 12.41 and 12.51 the Epstein zeta function is, up to a multiplicative 
factor of (27r) _2s , the spectral zeta function associated to the real torus RTa with A — Q*, where 
Q* is the dual lattice to Q. It can be shown that Z(s, Q) admits a meromorphic continuation and 
the functional equation 

Tr- s T{s)Z{s, Q~ l ) = (dct Q) 1 / 2 TT s - d ^ 2 T(d/2 - s)Z{d/2 - s, Q). 

The Iwasawa decomposition of Q asserts that Q can be uniquely expressed as 



Q 



1 o 

-*X Id-: 











with y G R+, x G R d_1 and Y a (d — 1) X (d — 1) positive definite matrix; here, Id-i denotes 
the (d — 1) X (d — 1) identity matrix and we think of x G R d_1 and y G R as row vectors. For 
m G R d_1 , one defines 

Q{m} = m ■ x + i\/yY(m), 

which is a complex number with Im(Q{m}) > unless m = 0. With all this, the generalization 
of the classical Kronecker limit formula for Z(s, Q) is the identity 



Z'{Q,Q) = -2-K^yZ{l/2,Y)-\og 



(2vr) 2 y J] t 1 



c 2-niQ{rn} 



(22) 



In [5] the authors study exp( — Z'(0, Q)) as a generalization of the classical Dcdckind eta function 
and prove an analogue of the Chowla-Selberg formula which relates exp(— Z'(0, Q)) to special 
values of higher order Gamma functions and special values of certain L functions. We refer 
the interested reader to [5] for additional results regarding the fascinating number theoretic and 
automorphic aspects of (|2"2"|) . 

2.8. Mellin transform and inversion. Let f x (t) = e~ t,>2 I x (t/2). Known transform identities 
include the formula 



Jo t y/irT\x + 1 - s) 



The identity (|23[) is valid for —x < Rc(s) < 1/2. The inverse Mellin formula gives 
/*(*) = ^~ [ Ms)t~ s ds, for - x < a < 1/2. 

27rl J (a) 



(23) 



(24) 



The Mellin transform of a product is computed via the convolution of two transforms, yielding 
the computations 

.dt 



fxfy(s) 



t 





1 

2tH 

1 
2nd 



—. I fy(z)t-*dz 



(a) 



(It 



for-y<cr<l/2 



dz, for —x < Rc(,s) - a < 1/2 



(a) 



fx{s - z)f y (z)dz. 
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Therefore, we conclude that 

fjy(s) = I e- t I x (t/2)I y (t/2)t 



tT ,. T , „dt 



i 

1 f T(s-z + x)T(l/2 + z-s)T(z + y)T(l/2-z) 



2ni J {a) ^T(x + 1 + z - s)^T(y + 1 - z) ^ 
1 f T(s-z + x)T(l/2 + z-s)T(z + y)T(l/2-z) 



2mJ {lj) ttT(x + 1 + z - s)T{y + 1- z) 

Using the bounds I x (t) = 0(t x ) as t — > and e~ t I x (t) = 0(£~ 1//2 ) as t — > oo, the Mellin transform 
fxfy(s) is defined for -x - y < Rc(s) < 1, so (J5SJ) holds for < cr < Re(s) < 1/2. If x ^ or 
y =/= 0, we can take s = and write 

_ t _ i r r(x-z)r(i/2 + z)r(z + j/)r(i/2-z) 



e- t [4(t/2)Jv(*/2)] T = ^- / p7 ! -i _i_ W TT ^ ( 26 ) 

i 2ni J(i/4) 7rr(a: + 1 + z)r(y + 1 - z) 

In section I6T21 below, we will extend (|26|) to the case x = y = 0. 

2.9. Miscellaneous results. For any e > and z £ C, Stirling's formula for the classical 
Gamma function is the asymptotic relation 

logT(z) = (z — 1/2) log 2 — z + (1/2) log(27r) + 0(1/ z) as z — > oo provided | arg(z)| < 7r - s. 

In particular, one has for fixed a and 6, the asymptotic relation 

| °j = z Q - & + O (z"-^ 1 ) as z ^ oo provided | arg(z) | < vr - e. (27) 

Further terms in the Stirling's formula can be computing (see, for example, |13j), which then 
would imply further terms in the asymptotic expansion (|27[) . 
Throughout our work, we will use the elementary identity 

oo 

/If 
{e^-e' wt )— for all w G C with Re(iu) > 0. (28) 

o 

To prove this relation, one simply observes that both sides of the proposed identity vanish when 
w = 1 and have first derivative equal to 1/w. 



3 Zeta functions and determinants for discrete tori 

In this section wc study the Gauss transform of the trace of the heat kernel associated to a discrete 
torus. We begin with Lemma |3~T1 which recalls the theta inversion formula associated to a general 
discrete torus. From the theta inversion formula, we define the Gauss transform (Lemma 13. 2[) . 
Using the group periodization representation of the heat kernel, we decompose the theta inversion 
formula into two summands: The identity term and the set of non-identity terms. We study the 
Gauss transform of these two summands separately, ultimately arriving at an identity (Theorem 
13. 6p which expresses the determinant of the Laplacian on the discrete torus in terms of 7-Bessel 
functions. 

3.1. Lemma. Let 6^(t) be the theta function associated to the d- dimensional discrete torus 
NZ d \Z d , defined by 

M*) = E e_A,t 

where 

{Aj} = {2d — 2 cos(27rmi/ni) — • • • — 2 cos(2irmd/rid) : < m, < n,, for each i = 1, 2, . . . , d}. 
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Then for all t > 0, we have the identity 

e N {t) = v{N) n ^ 2t in rkJ m, 

KeZ d l<j<d 

where K runs over d-tuples of integers (fci, . . . , kd). 

Proof. We refer to [TS] for the case d = 1. From there, one gets the general case using that the 
heat kernel on a product space is equal to the product of the heat kernels, since the Laplacian on 
the product space is defined to be the sum of the Laplacians from each factor space. □ 

3.2. Lemma. For all s G C with Rc(,s 2 ) > 0, we have 



oo 



2 2 \ =V(N)2s f e- s2t e- 2dt (I (2t)) d dt + 2s f e~ s2 * \o N (t) - V(N)e- 2dt T (2t) d - 1 



dt. 







Proof. By the definition of 8n, wc can write 

J2 e ~ Aj * = V ( N ) (e- 2t I (2t)) d + [0 N (t) - V(N)e- 2dt I (2t) d - 1]. 

Now simply multiply both sides of this identity by 2se~ s t and integrate with respect to t on 
(0, oo). Asymptotic behavior of the integrands as t — > and t —> oo easily imply that the resulting 
integrals are convergent for s G C provided Rc(s 2 ) > 0. □ 

3.3. Lemma The function 

f(s)= Y l°g(s 2 + A,) 
is uniquely characterized by the differential equation 

A^O ^ 3 

and the asymptotic relation 

/(s) = (F(iV)-l)-logs 2 + o(l) os s ^ oo. 



Proof. The differential equation characterizes /(s) up to an additive constant, which is uniquely 
determined by the stated asymptotic behavior. □ 

3.4. Proposition The function 

oo 

Us) =-[ (e- s2t e~ 2dt Io(2t) d - e~*) j 
o 

is uniquely characterized by the differential equation 

oo 

d s l d (s) = 2s J e- s2t e- 2dt I (2t) d dt 
o 

and the asymptotic relation 

Id{s) = log s 2 + o(l) as s — ► +oo. 
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Proof. For this, we write 

Us) = - / e-V 2 * ((I (2t)) d - l) j - I (e^e" 2 ^ e"*) | 



oo oo 

dt f / _2 t a 



o o 

By equation pg]). 

OO 

(e^e" 3 * - e" 4 ) - = - log(s 2 + 2d), 

= - I e-° 2t e- 2dt (I (2t) d - 1) j + log( S 2 + 2d), 
o 

From this last expression, the asymptotic behavior as s — > oo is immediate. □ 
3.5. Proposition The function 



so then we have 



n N (s) = - J (e- s2t [9 N {t) - V(N)e- 2dt I Q (2t) d - l] + e"*) 
o 

is uniquely characterized by the differential equation 

oo 

d s H N (s) = 2s { e- s ' 2f \e N {t) - V(N)e- 2dt T {2t) d - 1 a 



dt 

T 





and the asymptotic relation 

Wtv(s) = — log s + o(l) as s — > +oo. 

Proof. For this, we write 

oo oo 

?M*) = - | e" s2t (fjv(t) - V(N)e- 2dt I (2t) d ) j + J ( e - 2f - e"*) * 



oo 

= - / e~ sH (9 N (t) - V(N)e- 2dt I (2t) d ) j - log( s 2 ). 



(29) 



Clearly, the integral in (j2"9")l approaches zero as s approaches infinity, which completes the proof 
of the stated asymptotic relation. □ 

3.6. Theorem. For any s G C with Re(s 2 ) > 0, we have the relation 

J2 log (s 2 + A,-) = V(N)l d (s) + H N (s). 

A j7 t0 

Letting s — > 0, we have the identity 



where 



log [ [] a j ■ ) =v(N)i d (o) + n N (o) 



(e-^Io&f-e-^j 
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and 

oo 

Hat(O) = - J (8 N (t) - V(N)e- 2dt I (2t) d - 1 + e- f ) j. 
o 

Proof. We begin with the relation from Lemma l3.2l By substituting from the differential equations 
from Lemma l3.3[ Proposition 13.41 and Proposition 13.51 and then integrating, we get that 

J2 log (s 2 + A,) = V(N)J d (s) + H N {s) + C (30) 

for some constant C. We now use the asymptotic behavior relations as s — * oo from Lemma 13.31 
Proposition [3T4] and Proposition 13. 51 to show that C = 0. From the series expansion (fT3|) . we have 
that 

e~ 2dt h{2t) d - e~* = 0(t) as t -* 0. 

Lemma 14.11 gives that 

e- 2dt I {2t) d = 0(t- d ) asi^oo, 

so then the integrand in the definition of ld(0) is i 1 (0, oo) with respect to dt/t. Concerning the 
integrand in the definition of 7Yjv (0) we have that 

6 N {t) - V(N)e- 2dt {I {2t)) d = 0(t) as t -> 

and 

e"* - 1 = 0(t) as t 

so then the integrand in the definition of 7Yat(0) is in L 1 (0, 1) with respect to dt/t. Furthermore, 
we have 

0jsr(t) — 1 = 0(e~ ct ) as £ — > oo for some c > 

and 

y(7V)e- 2rft (/ (2t)) d - e~* = 0(t- d ) as i ^ oo, 

so then the integrand in the definition of Hn{0) is in oo) with respect to dt/t. With all this, 
we have that all functions in (|30[) are continuous and well-defined for s G R>0: so we simply need 
to evaluate at s = to complete the proof. 

□ 

4 Bounds and asymptotic formulas for Bessel functions 

In this section we prove bounds for individual /-Bessel functions. For technical reasons, it is 
necessary to separately establish bounds for Iq and for I x for x > 0. As will be seen in the next 
section, it is necessary to consider the function 

we" /„(„) (u t) 

for u > 1 and where n(u)/u — > x. The asymptotic behavior as u — > oo is given in Proposition 
14.7) uniform bounds are established in Lemma 14. II when x — 0, and in Lemma 14.41 and Lemma l4.6l 
when x > 0. Although we restrict our attention to integers x > 0, we recall that I x — I— x , so the 
results we prove apply for all x € N. 

4.1. Lemma. For any e < tt/2 and t > 0, we have the bounds 
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Proof. The positivity of Jo (t) follows immediately from the series expansion stated in section 
From the integral representation, choose any e £ (0, 7r/2) and write 



J («) = - / e-*( 1 - cos ("»d M +- / e -' (1 ~ cos(tl) W 



There exists c(e) > such that 1 — cos(u) > c(e)w 2 for u <E [0, e]. With this, for the first integral 
we have 

e e 



7T 



y / 4c(e)7rf 



For the second integral, we trivially have 



I /V* (1 - COs(tl)) dw < - E e -*(l-<=<»(e)) < 7r ~ g c -c(e)te 2 



Combining, we have 

< e-%(t) < , 1 + !^£ e - c (^ 2 . 
N /4c(e)7ri tt 

Since e < 7r/2 and u < e, we have the bound 

1 - cos(u) > u 2 /2 - u 4 /24 > (1/2 - e 2 /24)u 2 , 
so we may take c(e) = (1/2 — e 2 /24). Therefore, 

< e~%(t) < 1 + !L^£ e -(l/2- E 2 /24)e 2 ^ 

V(2-£ 2 /6)^ 7T 

Using elementary calculus, one shows that 

t l/2 e -ot < 1 / v ^ > 

so then 

< e-*/ (t) < C • r 1 / 2 where C = — 1 + ■ -== 1 (31) 

V(2-£ 2 /6V t v /(l-£ 2 /12)e 2 e 

□ 

4.2. Remark. Directly from Lemma I4TT1 we have that 

ue- uH I {u 2 t) <u-C- {uH)- 1 ' 2 = C ■ t- 1 ' 2 , 

so wc indeed have established the uniform upper bound as claimed. Also, for our purposes, it 
is not necessary to optimize pip through a judicious choice of e. Numerically, one can show 
that by taking e = 1.5, which is allowed since we only required that e < tt/2, we have that 
C = 0.676991 .... One point we will use later (see section l7TTj) is that C < l/y/2 < 1. We will use 
the numerical verification of this bound and omit the theoretical proof from our analysis, noting 
that the estimate for C indeed can be proved from (f3"Tj) . 



4.3. Lemma. For fixed x > 1 and t > 0, consider the function 

t 



h x (t) = v/(x 2 +t 2 ) -t + xlo? 
Then we have the bound 



exp(h x (t)) < 

By continuity, the inequality also holds when t = 



t 



x + ^{x 2 + 1 2 ) ^ 

x/2 



t + X 



Proof. To begin, observe that 

„2 



h x(t) = 7= = +x\og 



t+^{ X 2 +t 2 ) \x + yj{x 2 +t 2 ) J ' 

which comes from the definition of /i x (t) ancl by writing 



Therefore, we have that 



h x {t) = , + 2: log 





ft 2 ) 












-t 2 ) 








-t 2 ) 


x 2 





.r 2 



f + v/^+t 2 ) 



^+ ^(x 2 +t 2 ) ) 
x log ^- {^J X 2 + t 2 — x^j ^ 



a: log J 1 + — -- 



t + v/(x 2 + t 2 ) °VW *v * 



We now employ the change of variables 



which is equivalent to the relation sinh(u) = x/t. Using the elementary identities 
1 + (sinhtt) = (cosh it) 2 , sinhw = 2sinh(u/2) cosh(«/2) 

and 

1 + coshu = 2(cosh(u/2)) 2 , coshu — sinhw = e~ u 
to arrive at the expression 

h x {t) = x (tanh(u/2) — u) . 
Trivially, since it > we have that tanh(u/2) < u/2, so then, for x > 0, 

h x (t) < -x ■ u/2 = -(a/2) log (J (l + ^ + ^ < -(x/2) log (l + 
With all this, we have that 

>(M*)) < (! + f) 



-as/2 / t ^ X/2 



exp( 



x + 1 



which completes the proof of the lemma. 

4.4. Corollary: For any t > and integer x > 0, we ftave 

^-'«" s (^r=( i+ Tr /a 
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Proof. We begin by considering x > 2, so then the analysis from Lemma T4. 31 applies. Indeed, we 
use the trivial estimates 

1 =e l/(2V(* 2 +t 2 )) < < 1 



V2tt V2tt 

and 



(X 2 +< 2 ) 1 /4 

Therefore, using the notation of Lemma 14. 3[ the bound (fT5|) becomes 

a/2 



Vt-e-H x {t) <exp(M*)) < ( 



which proves the claim, again provided that a: > 2. If x = 0, the claim follows from Lemma |4~T1 as 
well as Remark 14.21 which shows that the constant C in Lemma |4~T1 satisfies C < 1. It remains to 
consider the case when jc = 1. 
The series representation (fT3| of (f ) gives 



For i < 1, i/2 < 1/^(4+1), so then 

1/2 



for t < 1. 



From ([T4|). we have that Ji(i) < -fo(i) since cos(x6 | ) < 1. With this, Lemma \4. II implies 

e~ l h{t) < e-*I (t) < C/Vt (33) 
In Remark 14.21 it was argued that C < 1/V2- With this, we have for t > 1 the inequalities 



c/Vt < 1/V2I < W(t + i), 

so then 

e"*/i(i) < 1) fort>l. (34) 

Combining (|33|) and (|34|) . we have the claimed assertion for x = 1, which completes the proof for 
all integers x > 0. □ 



4.5. Remark. To be precise, the bound in Corollary 14.41 in the case x = 1 is not needed in this 
article. We included the statement and proof for the sake of completeness. 

4.6. Lemma. Fix t > and non-negative integers x and hq. Then for all n > no, we have the 
uniform bound 

/ . \ nox/2 / \ -nox/2 

0<^.e-" 2 U 2 ()< M^- =( 1 + A) 

\x + n tj \ n tj 



Proof. The positive lower bound is obvious from the series definition of the /-Bessel function, so 
we focus on proving the upper bound. From Corollary 14.41 we have 

v/2 



^J(nH)-z- nh In X {nH) < ( nH 

\nx + n t 



Elementary algebra yields 



nH \~_ 1 x_ = 1 , *7(2t) 



nx + n 2 t J nt (nx/2) 
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so then we have 

/7-5-r nH , as ( riH \ nx/2 ( x 2 /(2t)\' nx/2 



(35) 



For any constant c > 0, consider the function 

, y 



g{y) = [ i + c 



We claim that g(y) is monotone increasing in y. Indeed, using logarithmic differentiation, we have 
that 

-logfl + ^+^^V~ = logfl + -V^=log YC + 2/ 



g(y) \ y) l + c/y y 2 \ y) c + y \ y J c + y 

= -logP^V — = -log(l-^-V— ■ 
\c + yj c + y \ c + yj c + y 

Clearly, if < u < 1, the function — log(l — u) — u is positive. Since g(y) > 0, we conclude 
that g'(y) > 0, which shows that g(y) is monotone increasing in y. By taking c = x 2 /(2t) and 
y = nx/2, we have, for n > no, the inequality 

nox/2 / 2 //o-A\ 



^my ^" ( x 2 /{2t) \ 

+ (nox/2) J " V (W2)/ 



which gives 

1 + x 2 /(2t) \ - nx/2 < / z 2 /(2t) \ ~" 03;/2 _ / -" oa:/2 _ / npt N " oa;/2 

(nx/2) J ~ \ (n x/2)J \ n t J \n t + x 



(36) 



By substituting (|36|) in (l35|) . the lemma is proved. □ 

4.7. Proposition. Let n(u) be a sequence of positive integers parameterized by u £ Z + suc/i £/ia£ 
n(u)/u — > a > as oo. TTien /or any t > and non-negative integer k > ; we /iai>e 

lim n(u) • e- 2 " 2t /„ (u)fc (2u 2 t) = • e ^ 2 '^. 



Proof. We work with the integral expression for the /-Besscl function given in section I2T21 namely 
that for any integer x, we have 

I x (t) =1 J e tcos(e) cos(0x)d8. 
o 

To begin, assume k > 0. If we let y = n(u)k9, then we can write 

7T 

n(n)e- 2 " 2 */ nHfc (2n 2 = ^ J e -^{i-co S {e)) cos ( 9n ( u)k)d0 

o 

n(u)fc-7r 

1 



kit 



e -2u"t(l-co B (»/(n(u)fc))) cos ( 1/ j dj/- 



For all v G [0, 7r], one can easily show that 



1 7r 2 1 w 2 1 — cos(v) 
2~24-2^24- ~ 2 ' 
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Let c = 1/2 — 7r 2 /24, which, numerically, can be shown to satisfy c > 0. Setting v = y/(n(u)x), 
we get for y € [0, ux7r] the uniform bound 

1 - cos(y/{n(u)k)) / 



1/u 2 \n(u)k 

In addition, observe that 

1_- cos(y/(n(u)k)) _l ( y \ 



l/u 2 2 Vafc- 

Choose Mo so that for u > uo, we have a/2 < n(u)/u < 2a. Using elementary bounds, we have 
for any u > uq the inequalities 



n(u)k-K 

/net 

Therefore, by the Lebesgue Dominated Convergence Theorem, we have 



/ s -2u 2 t r /r, 2j.\ ^ 1 / -cy 2 t/(2ak) 2 j ^ 1 / -cy 2 t/(2ak) 2 i 

A u ) e I n (u)k{*u t) <— I e " IK ' dy < — I e y n > dy 



hm n(«) • e- 2 " 2t / nHfe (2u 2 i) = J- ) e^l^f C0S (y)d y = -J= • e -(" fe ) 2 /( 4t ). 



If a; = 0, the proof follows a similar pattern using instead the substitution y = 9u; for the sake of 
brevity, we omit the details. □ 

4.8. Remark. Proposition ^. 71 was proved by Athreya p] Theorem 2] using a certain local central 
limit theorem. For convenience we included a quick self-contained proof. 



5 Asymptotic behavior of spectral determinants 

In this section we use the bounds and asymptotic relations from the previous section to prove 
our main theorem, namely the asymptotic behavior of the spectral determinant associated to a 
sequence of degenerating discrete tori. We begin by proving that the associated family of traces of 
heat kernels converges pointwise through degeneration (Proposition 15. 3| ) and then prove uniform 
bounds for heat traces for long time (Proposition [5T5]) and small time fProposition l5.6l and Lemma 
15. 7|) . After these results, we will analyze the expression for Ti.^(0) from Lemma l3.5[ namely 

OO 

H N (0) = - J (e N (t) - V(N)e- 2dt I (2t) d - 1 + e-*) j. (37) 
o 

Ultimately we compare the limiting value of (|37p through degeneration with the expression for the 
spectral determinant on the real torus, which is stated in (|2ip in section \2~M These computations 
are given in the proof of Theorem 15.81 which is the main result of this section. 

5.1. Degenerating sequences of discrete tori. We consider d-tuples of integers N(u) 
parametrized by a positive integer u in such a way that 

-N(u) = ~(m, ■ • • , n d ) — > (a 1: ■ ■ ■ , a d ) as u — > oo. (38) 
u u 

Let A be the diagonal matrix with the ay's on the diagonal. Recall the notation V(N) = rii ■ ■ ■ rid 
and V(A) = a.\ ■ ■ - ad- From (|38|) , V(N{u))/u d — * V(A) when u — > oo. For the sake of brevity, 
and, for this section, we assume that V(A) > 0, meaning the limiting real torus AZ d \R d has 
dimension d. Using the elementary change of variables 1 1— > u 2 t, we write (|37j) as 

oo 

n N(u) (0) = -J ^ JV(u) (u 2 i)-U(iV( U ))( e - 2 " 2t /o(2u 2 t)) d -l + e -" 2t ) j, (39) 
o 
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which is the form that we will study. The results of this section are designed to determine the 
asymptotic behavior of (|39[) as u — ► oo. 



5.2. Proposition. For each fixed t > 0, we have the pointwise limit 

#jv(u)(w 2 i) — > Oa(£) asw^ oo. 

7n words, the rescaled theta functions on the discrete tori approach the theta function on the 
limiting real torus AZ d \R d . 

Proof. We begin by writing the theta function using its expansion involving 7-Bcsscl functions, 
namely 

oo d 

6 N{u) (uh)= Y[n ] (u)e- 2u2t I nj{u)k] (2u 2 t). (40) 

From Proposition 14.71 we have, for any t > 0, the pointwise limit 

d d 1 

j=i j=i V I ) 

Since 

oo d 
fcl,". ,k d = -oo j = l V V ^ 

the result will follow if we can interchange the limit in u with the infinite sum in (|40[) . 

Choose uo such that for each j = 1, • • • , d, otj/2 < rij(u)/u < 2a j for all u > uq. We can re- write 
(|40p into a sum of d + 1 subseries, each determined by the number of fc/s which are equal to 
zero. Fix some uo > sufficiently large. Then from Lemma |4 . 1 1 and Lemma l4.61 we have that the 
subseries consisting of terms with exactly r of the fcj's equal to zero is bounded from above for 
any u > uq by 



oo ' / h \ ~' u o a jkj /2 oo /■ 

2 r C r {2t) -r,2y {A) £ Tj(l + gj < 2^(2*)-^^) £ ft rj 



A:i j--. .k r — 1 j — 1 v k\ — 1 j — 1 

(41) 

where 

-tlOQj/2 

< 1. (42) 



4u t 



Obviously, 

oo r r / 

£ IIv 1I(t^- 

fci,-..,fc r =l3=X 3=1 V J 

Therefore, the series in (|40p is uniformly convergent for fixed t, so we can interchange the limit in 
u and the summation in (|40p . which completes the proof. □ 



5.3. Lemma. Let 

oo 

#ab S (0 = 25>- CJ ' 2t 
3=1 

with c = 4ir 2 (l — 7r 2 /24) 2 . Assume uq is such that a k /2 < n k /u < 2a k for all k = 1, • • • ,d and 
u > uq. Then for any t > and u > uq, we have the bound 



N {u){u 2 t) < n (l + e" 4 ""* + 9 &hs {a k t/2) 



fe=i 
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e -4(sm{iTj/n k (u))yt 



Proof. Recall from [T5] that 

n k (u)/2-l 

Q nk(u) (t) = l + e- At + 2 J2 

i=i 

if n k (u) is even, and 

(n fe («)-l)/2 

0n k (u)(t) = 1 + 2 e - 4 ( sin (^'/^(«))) 2 * 
3=1 

if n k (u) is odd. In cither case, we will use the elementary bound 

smx>x — x 3 /6 for x £ [0, 7r/2]. 

To prove this inequality, one considers the Taylor series expansion with error at x = for sin(x), 
together with the observation that the error term is positive. Taking x = ?rj /n k (u), we then have 
that 

n k (u)sm(irj/n k (u)) > - (irj) 2 / (6n k (u) 2 )). 
We may assume that j < n k (u)/2, so then j/n k (u) < 1/2, hence 

n k (u)sm(irj/n k (u)) > 7T?'(1 - 7r 2 /24). 

It is important to use the trivial bound 1 — 7r 2 /24 > 0. With all this, we arrive at the inequality 

n/2 oo 
e nj(u) (u 2 t) < 1 + e- iu2t + 2^2 e" ca ^ 2 */ 2 < 1 + e" 4 ^ 1 + 2^2 e" ca ^' 2 * /2 

3 = 1 3 = 1 

with c = 47r 2 (l — 7r 2 /24) 2 > 0. Therefore, we have shown that for any u > uq, we have the bound 

e n {uh) < i + e - 4u °' + # abs (i). 

Since 1 

d 

3=1 

the proof of the lemma is complete. □ 
5.4. Proposition. With notation as above, we have that 

(6 N(u) {u 2 t) - V(N(u))(e~ 2u2t I (2u 2 t)) d - 1 + e-'" 2t ) 



dt 

T 



(e A {t) - 1) j - ^V{A){Air)- d ' 2 + o(l) as u -> oo. 



Proof. Write 



^ (M) (u 2 t) - y(7V( U ))( e - 2 '" 2 %(2 U 2 t)) rf - 1 + e"" 2 ') j 



OO OO 

di /" -2u 2 t T ,o„.2+\\d dt , [ „~u 2 t dt 



(43) 



(fljv W (« 2 «) V(N(u))( e -^ %(2uH)) d j + j e 



21 



and we consider the three integrals separately. Trivially, we have 



oo 



-„ 2 t dt 

e > as u — ► oo. 

t 



Next, we claim that 

oo 

, d dt „, /\. , j/,di 2, 



y(JV) / (e- 2tl2t / (2u 2 t)) y -> V(A) j (Antr^j = ^V(A){4m)- d / 2 as u - oo. (44) 
l l 

Indeed, the pointwise convergence of the integrands in (|4"4"|) is proved in Proposition 14.71 and the 
uniform bound from Lemma 14.11 which is integrable on (1, oo) with respect to the measure dt/t, 
allows one to apply the Lebcsguc Dominated Convergence Theorem to prove the claim. It remains 
to show that 



(e N[u) {uh)-i) d l^ J (Q A (t)-i)j 
i 



as u — > oo. 



The pointwise convergence of the integrands is proved in Proposition 15.21 and Lemma 15.31 estab- 
lishes a uniform, integrable upper bound so that, again, we may apply the Lebesguc Dominated 
Convergence Theorem, thus completing the proof. □ 

5.5. Proposition. With notation as above, we have that 
l l 

(e N{u) (u 2 t) - V(N)(e- 2u2t I (2u 2 t)) d ^ j -> f (e A (t)-V(A)(A7rt)- d / 2 )j 



o 

as u — ► oo. 



Proof. For fixed t < 1, we have the pointwise convergence 

e N {u){uH) - V(N)(e- 2u2t I (2u 2 t)) d Q A {t) - V(A){AKty d ' 2 as u 



directly from Proposition [477] and Proposition ^. 21 It remains to show that we have uniform, inte- 
grable bounds for the integrand so then the proposition will follow from the Lebcsguc Dominated 
Convergence Theorem. For this, we write 

e N{u) {uH)-V{N{u)){e- 2uH I {2uH)) d = V{A)Y, J[ ue- 2u2 %. (u)kj (2u 2 t). (45) 

The bounds ([4Tj) and (|42j) apply. Since each rij(u) tends to infinity, one can choose u$ so that 
for each j = 1, • • • , d and u > uq, we have the inequality rij(u)kj > d + 2, which implies that the 
upper bound in (j4"2"|) is integrable on (0, 1) with respect to dt/t. □ 

5.6. Lemma. For u G R, we have the asymptotic formula 
l 

J (e- uH - l) j = T'(l) - log(u 2 ) + o(l) as u -> oo. 
o 

Proof. Choose e £ (0, 1). By employing integration by parts, we can write 

11 l 

- u2t - l) j = J (e-" 2t - l) dlog(t) = log(t) • (e"" 2 * - l) ||_ + u 2 J \og(t)e- u2t dt 

S € 

1 

= log(e) • (l - e-" 2£ ) + u 2 J log(t)e- uH dt. 
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For fixed u, we have that 1 — e " e = 0(e) as e approaches zero, so then 

1 l 

e -« 2 * _ l) ^ = u 2 f log(t)e- u2t dt. 



o 



If we let v = u 2 t, then dv = u 2 dt, hence 

1 u 2 

dt 



(e"" 2 ' - l) — = / \og{v I u 2 )e' v dv = / \og(v)e- v dv - \og(u 2 ) / e~ v dv 



By definition 

oo 

, dv 



r( s ) = / e-v 

which implies that 



r'(l) = j e~"log(v)dv = j e~"\og(v)dv + / e" u \og(v)dv, 





thus yielding the relation 



1 oo 

dt 



l) — = r'(l) - / e~ v \og{v)dv - \og(u 2 ) / e~ v dv 



o 



(46) 



= r'(l) - log(u 2 ) + log(M 2 ) e - M - J log(v)e- v dv. 

U 2 

For x > 1, we have the trivial inequalities < logx < 2e x / 2 , so then 

oo oo 

< / log(v)e- v dv < 2 / e-" /2 dv = 4e-" 2/2 = o(l) as u -> oo. (47) 



Trivially, log(u 2 )e " = o(I) as u — > oo. By substituting this estimate and (|4"7]) into (14"6")) . the 
proof of the lemma is complete. □ 

5.7. Corollary. With notation as above, we have 



l 



(o N{u) (u 2 t) - V(N) (e- 2 ^I {2u 2 t)) d 1 + e-^ j 

[0 A (t) - V(A)(4irty d / 2 ) j + T'(l) - log( U 2 ) + o(l) 



dt 

as u — ► oo. 



Proof. Simply combine Proposition 15 . 51 and Lemma 15.61 □ 
5.8. Theorem. With notation as above, we have 

Ww(«)(0) = logu 2 - Ca(°) + as u -> °°- 



23 



Equivalently, we have 



log [ J] A 3 ; ] - V(N(u))l d (0) + logu 2 - Ca(0) + o(l) as u oo 



where 



U0) = -[ (e- 2dt I {2t) d -e~t)*. 
o 

Proof. To begin, combine Proposition 15.41 and Corollary 15.71 to get the asymptotic expansion of 
V(N(u))TCn( u <) (0) out to o(l). Then combine the terms in the expansion with the expression in 
(l2"Tj) to complete the proof of the first assertion. One then substitutes the first assertion into the 
identity stated in Theorem 13.61 to prove the second relation. □ 



5.9. Remark. Using the notation of spectral determinants, one can re- write the main result in 
Theorem 15.81 as stating the asymptotic relation 

logdet*A £)TiA r( tl ) = V(N(u))ld{0) + logu 2 + logdet*A flT: A + o(l) as u — > oo. 

In words, the asymptotic expansion of the log-determinant of the Laplacian on the discrete torus 
DTn^ has a lead term 

oo 

V(N(u))X d (0) + logu 2 = -V(N(u)) I (er 2dt (I (2t)) d - e"*) j + logu 2 



V(N(u)) ( log 2d- / (e- dt (I Q (t) d -l) j\+l 0i 



a 2 



with a constant term equal to the log-determinant of the Laplacian on the real torus AZ d \R d . 

5.10. Remark. As previously discussed, there are two zeta functions being considered: The 
spectral zeta function (a associated to the real torus AZ d \R, d and the Epstein zeta function 
Z(-, Q) associated to the positive definite d x d matrix Q. Classical and elementary mathematical 
considerations show that Z(-, Q) = (a where Q is the form associated to A* , the dual lattice of A. 
This remark needs to be kept in mind when reading the main result in Theorem 15.81 as restated 
in Remark 15.91 



6 Example: the cases d = 1 and d = 2 

As stated in the introduction, the lead term in the asymptotic of logdet* A. DT N ^ for d = 2 was 
first computed by Kasteleyn |16| . Duplantier and David [S] expressed the next order term in terms 
of the Dedekind ry function. In this section we rederive these results by the methods introduced 
in this paper and go further by indicating how to obtain the complete asymptotic expansion of 
logdet* A.]jt,n(u) when d = 2. We begin by making explicit the case d = 1. 

6.1. The case d = 1. Known evaluations of the Riemann zeta function (q at s = yield the 
relations Cq(0) = -1/2 and Cq(0) = -(1/2) log(2?r). Let N(u) = n with n/u -> a > 0. With all 
this, we can write 

logn 2 =logw 2 +41og(27r/a)CQ(0)-4Cq(0) + o(l) as u -> oo. (48) 

We claim that (|4"5)) is equivalent to the Main Theorem when d = 1. First, it is obvious that the 
number of spanning trees in nZ\Z is n and, hence, the determinant of the combinatorial Laplacian 
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is n 2 , so the left-hand-side of ([5]) is log n 2 . Regarding the right-hand-side, it is possible to directly 
verify through the Mellin transform that 



oo 

J e- 2t (I (2t)~l)^ = log2, 



which implies, in our notation, that Xi(0) = 0. Hence, the lead term in the asymptotic expansion 
([5|) when d = 1 is zero. The continuous Laplacian on nZ\R has the set of eigenvalues given by 
{(2ir) 2 (n / a) 2 } for n <G Z, so the spectral zeta function is 

Us) = 2(2Va)- 2s C Q (2s), 



so 

C(0) = -logdet*A flT , Q = -41og(27r/a)C Q (0) +4^(0), 
which agrees with (|48p and confirms, by direct computation, the Main Theorem when d = 1. 

6.2. The lead term when d = 2. We will evaluate the expression 

oo 

I d (0) = I log(2d) - / {e- dt (I (t) d - 1) j 



when d = 2. Recall from section (|2.8|) the Mellin inversion formula, namely equation (|25|) . Let us 
set x = y = and move the contour of integration to the line Re(er) = —1/4. The integrand has 
a pole at z = 0, which results in the expression 



~ i r r( s -z)r(i/2 + z- 5 )r(z)r(i/2-z) j r( s )r(i/2 - s)r(i/2) 
/o/o(s) = ^— / .tvt , : ^ — ^ dz + 



2m J nT(l + z - s)T(l - z) 7iT(l - s)T(l) 

Re(z) = (-l/4) 

so then, using the definition of T(s) as the Mellin transform of e - *, we have 

/*°° 

/ / ( s )-r( s )= / e -'[/ (i/2) 2 -ir^ 

Jo ' 

1 f T(s-z)T(l/2 + z-8)T(z)T(l/2-z) 



(49) 



2tti 7 7rr(l + z - s)r(l - z) 

Re(z) = (-l/4) 

r( 8 )r(i/2- s ) 

+ — i=rn ^ r s ' 

^^(1 - s) 

By Taylor's theorem, we can write 

hm - r (s) ) = -rW2) + ^(D 

s-o ^ 0FT(1 - a) V 7 
where we have used that T(l/2) = y/% and T(l) = 1. One now uses the relation 

T(z)T(z + 1/2) = 2 1 - 2z y ^T{2z) 
and, by computing the derivative at z = 1/2, one can show that 

v/tt 

With these preliminary results, we have shown, by taking s — ► in (|50p . that 

, Ir , )W2 ,7* 1 /" r(-z)r(i/2 + z)r(^)r(i/2-z)^ 



(50) 



e -*[7 (t/2) 2 - 1]- =— / v ; ')' ■ '-dz + log4 

t 2m J 7rl (1 + z)l (1 — z) 

Re(a)=(-l/4) 

1 f 1 (51) 



iog4 -^ y 

Re(»)=(-l/4) 
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where the last equality of integrals comes from employing the identities 
r(a;)r(l — x) = 7r/ sin7ra; and xT(x) = T(x + 1). 

We now evaluate the last integral in (|51[) by moving the line of integration toward Rc(z) = — oo. 
The integrand has poles z = —n — 1/2, for n = 0, 1, 2... with residues (— l) n (n + l/2)~ 2 /7r. Thus 

,dt 1 ^ 4G 



**,,■,■> , .dt 1 -1)" 4G , N 

where G is the Catalan constant 



With all this, we have shown that 



(2n+l) 2 
4G 

X 2 (0) = — . (53) 

The relation in (|53[) . which is the main asymptotic term from Theorem 15.81 in the case d = 2, 
agrees with the computations from [9], specifically equation (3.18) on page 349; see also equation 
(A.21) on page 427 of [9]. 

6.3. Secondary terms when d=-2. From Theorem 15.81 we have that 

H N (0) = logu 2 + logdet*A i?TidiA + o(l) as u -► (54) 

for a degenerating sequence of discrete tori. As before, consider a sequence of integer vectors N(u) 
for u G Z + and assume 

—N(u) — — (m(it), 712(1*)) — > (ai, 0:2) as u — > 00. 
u u 

Proposition 14. 71 evaluates the heat kernel on (a±Z x a 2 Z)\R 2 where the Laplacian is 

d 2 d 2 



A 



R 2 



dx\ dx 



where x\ and X2 are the usual global coordinates on R 2 . The action of Ar2 on the space of smooth 
functions on [a,\L x a 2 Z)\R 2 has the set of eigenvalues given by {(27r) 2 ((n/ai) 2 + (m/a 2 ) 2 )} for 
n, m G Z, so then its associated spectral zeta function is 

(„,4#(0,0) (K«l) 2 + W«2) 2 ) S 

Let us write 

Ct(s) = (aia 2 ) s £ , (z,s) 

where 

E(z,s) = (2tt)- 2s V 1 with z = i(a 2 /ai) = iy. 

z — ' \nz + Tn\ 

The function E(z, s) admits a mcrmorphic continuation to all s G C with expansion 
E(z,s) = -1 - slog{y\r](z)\ 4 ) + 0(s) as s -► 

where 2 = i(a 2 /ai) = iy and ry(z) is the classical Dcdekind eta function 

00 

rj( z ) = e 2mz JJ (l - e 27rmz ) for any z G C with Im(z) > 0. 

n=0 
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With all this, we then have that 



logdet* A RTtA = -ttCa = log(aia 2 ) + log(y|?7(z)| 4 ). 

OS s=0 

Therefore, 

H N (0) = logu 2 + Iog(aia a ) + log(y ^(2;) | 4 ) = Iog(m(«)na(u)) + log(y|77(z)| 4 ). (55) 

The expression (|55p agrees with the second-order term given in equation (A. 21) on page 427 of 
[9], where those authors used the notation M = ni(u) and N = n,2{u). 



6.4. Error terms when d = 2. The asymptotic expansion in Theorem 15.81 contains an error 
term of o(l). We present here a technique which improves the error term, ultimately obtaining 
an expansion with error term of arbitrarily small polynomial order as stated in the introduction. 
Additionally, the computations we give here give a second proof of the Main Theorem in the case 
d = 2. 

Recall the Mellin inversion formula from section (|2.8|) . namely 



ts dt_l f r(a -z + 3QT(| + z- s)T(z + y)T{\ - z) 



IJt/2)lJt/2)e-H s — = — / — —. ' 2 ^-Af Ldz. (56) 

w ' vy 1 ' t 2m J kT{x + 1 + z - s)T(y + 1 - z) V ' 

Re(z)=er 



For fixed a and b, Stirling's formula yields 

+ O^ " 6 - 1 ) (57) 



r ( z + «) _ a-b , 



T(z + b) 

as z tends to infinity anywhere in the half-plane Rc(z) > S > 0. Substituting ((57|) into (|56|) . we 
have that the lead term in the asymptotic expansion of (|56[) is 



Ro(z) 



J x s-*-l-z+s yZ -l + z . I r(z + 1 _ s)r( l _ z)dz 

(58) 



x2S i / (?fr ( , + ±-.)r(i-,)* 



7r • xy 2ni 

Kc(z)—fT 

Let us analyze ([58| by moving the contour of integration toward Re(z) = 00. The integrand has 
poles at z = 1/2 + n with residues equal to (— l) n /n\, so then (|5"5|) is equal to 

X 2s ^ /y\2„+l (-1)' 



E(*r'^»-. +1 ). o») 

Re-write (J5HJ) as 

y 2s-2 ~ /yx 2«-2s+2 (_l)r 



s= , ^r(„-« + i) 

n— 

and use the elementary formula 



T(w)a- W = I e~ at t w — (60) 



with w — n — s + 1 and a = x 2 jy 2 to get that 

00 

y 2s-2 ^ f^_ (xVv ^ n _ s+1 (-1)" dt 

n\ t ' 
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Observing that 

(-1)* 



-* = £** 



e 

we arrive at the relation 



S=tLl.f e -^^>e-H-^-. (61) 

IT J t 



The integral can be evaluated using ([60]) with w = 1 — s and a = (x 2 jy 2 + 1), yielding 

S = V— . T(l - s) (x 2 /y 2 + l) 3 - 1 = -T(l s) (x 2 + y 2 ) 3 - 1 . (62) 
If we set s = —w with Re(iy) > 0, we then have that 

oo 

/ I x (t/2)I y (t/2)e-H- w ^- = - ■ tg) + bwer order terms. (63) 

J in [x 1 + y 2 ) 

Finally, we set x = nj and y = mk with fixed j and A: and positive integers n and to and sum over 
all n and to, thus proving that 

oo 

J2 [ I nj m)I mk (t/2)e-H- w ^ = E^±^l - - j + lower order terms. (64) 

TZ.J 1 n \\ni) 2 + (mk) z ) 



a . in 



The series in (|6"4"|) can be related to the Eisenstein series Ea(z, s) with A = (x,y) and z — im/n. 
One then can use the classical Kronccker's limit formula to evaluate the asymptotic behavior of 
(|64p as w approaches zero. In all, the above computations provide another proof of the Main 
Theorem when d = 2. 

Stirling's formula can be used to compute further terms in the asymptotic expansion in (|57D . 
Specifically, for any positive integers N and M, we write 

TV T !Z Z) ^ Z+ J \ --(-)" = E ^Ws-^-V-" 1 -* + O (*— V*" 1 "") 
1 (# + 1 + z)l (y + 1 — z) xy \xJ v ' 

(65) 

as x,y — ► oo, where P^fc denotes a polynomial in z, and the sum is over all integers j and fc such 
that < j < N and < k < M provided (j, k) ^ (0, 0). By comparing with (|56p . we see that the 
lower order terms in (|64[) can be determined by studying 



/ ^W(f) 2Z r(z + i)r(i-z)dz. 



Re(z)=o- 

Moving the contour of integration, we get that 

2n+l 



I OO 



n=0 

Since Pj t k{z) is a polynomial in 2, we have that 

^ + i/»er-p«(il)(sr 

so then 

\ 00 
'y d \ f xy 



(66) 
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Let us define 

B 3,k(u) = ^ S jtk {niu,n 2 u). (67) 

(ni,n 2 )^(0,0) 

Elementary bounds show that 

Bj,k(u) = O (u~ j ~ k ~ 2 ) when u — > oo, provided (J, k) jt (0, 0) 
since the series in (|57|) is convergent for (j,k) ^ (0,0). Finally, if we set 

F„(u)= B i,k(u>), (68) 

j-\-k— n 



we have that the error term o(l) in Theorem 15.81 can be improved to 



K 



^ F„(u) + O (u K 3 ) as u — > oo for any integer A" > 0. 



n=l 



To summarize, the coefficients (|68f can be explicitly determined from the evaluation (|66|) which 
utilizes the precise evaluation of Stirling's formula as stated in (f6"5"l) . The coefficients (|55| can be 
made explicit by evaluation the polynomials Pj.k(z) as defined in We will not pursue this 

analysis further here, but instead leave the details for consideration elsewhere. 

7 Additional considerations 

In this section we present a number of computations and remarks discussing further aspects of the 
preceeding analysis. 

7.1. A precise evaluation of the lead term. Let 



oo oo 

a d = log(2d) - J e- dt (I (t) d - 1)* = log 2 - / {e- dt I (t) d - e" 4 ) 



dt 

T' 



which is obtained by replacing 2i by i in the integral definition of T d (Q). With this, 

oo 

a d+ i-a d = J e - dt (/ (i) d (l~ e -%(t))^. 
o 

Since Ia(t) < e 4 , we have that a d +i > ad- Let us write 

oo 

e-%(t) = l + ^M fe , 
fc=i 

which implies that 



oo 

^6„y e- dt (I (t)) d t k dt. 

From [15], we have the identity 



k=l 





Therefore, 



oo 

= / e-"*I (i)dt = — -L=. 



OO 

e- d *(7 (t)) d t fe ^ - (-9 u ) fc F* d (d), 



n 
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where F* d denotes the d-fold convolution of F with itself, evaluated at d. Therefore, we have 

00 

a d+ i =a d + Y J bk -(-d u ) k F* d (d). 



k=l 



In section [6721 we have shown that 02 = 4G/-7T, where G is the Catalan constant. In addition, one 
can write 



e-*/ (t) = - /e-* sin2(e/2) ^, 



so the coefficients bk are explicitly computable. With all this, we have presented another series 
expansion for ad which could be used to numerically evaluate ad- 

7.2. Numerical estimation of the lead term. This subsection should be compared with the 
discussion in Felker- Lyons [llj . By focusing on the lead term in Theorem 15. 8[ we have that 



00 

^ £ logA j (u) = log2d-Je- dt (I (t) d -l)j + o(l 



V(N( i 



) as u — ► 00. (69) 



As stated in the introduction, the explicit form of the eigenvalues are such that we recognize the 
left-hand-side of (f6"9"| as a Riemann sum, so then (f6"9")l implies the identity 

00 

J log (2d- 2cos(27ra;i) 2 cos(27rx d )) dx x ■ ■ ■ dx d = log 2d- J e~ dt (I (t) d - 1) j. (70) 

Z d \R d 

Recall that 

h{t) = 1 + a nt 2n where a n = 

n=l "•' 

Let us write 

00 

(I (t)) d = l + ^a„^ 2 ". 

n=l 

The coefficients a„^ easily can be written in terms of a n and multinomial coefficients, and (1701) 
can be evaluated as 



/oo 
log (2d — 2 cos(27ra;i) — • • • — 2 cos(27ra;d)) dxi • • • dxd = log 2d — 



a„ i(J r(2n) 



Z d \R d 



The numerical evaluation of the left-hand-side of (fTTTj) is difficult for large d, both because of the 
slow growth of the integrand and because the domain of integration involves the d-fold product of 
the unit interval. By comparison, note that the right-hand-side of (|T0[) allows for rapid numerical 
evaluation which improves as d gets larger. Indeed, from Lemma |4. II we have 



OO 00 00 j m t 1 

-rdj.-d/2+l dt 



< 
/(I 



Je- dt (I (t) d - l)j <C d J t-^ 2 j + i |e- dt dt = 2Ct ° ^ +6 ' ° . (71) 



This bound obviously gets closer to zero as to grows for any d. In addition, the bound decays 
exponentially in d provided Ct < 1 which, from the estimate that C < l/\/2, implies that we 
need to > 1/2. By the Mean Value Theorem, we have, for any z and a, the estimate 

(x + a) d -x d < ad{x + a) d - 1 . 

For any N, to be chosen later, we take 

N-l 



= S~] a nt 2n and a = I (t) - x < I (t) 



.r 

n=0 
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to arrive at the bound 



J e- dt (i (ty - 1)7 < / e ~ dt HE a " t2 "J ~ 1 )J +d I et ( e ' tJ «•(*)) ,, ' 1 (E ^* 2n J * 

Choose and fixed to € (1/2, 1). For any e > 0, there is an iV such that 

00 

Y,a n t 2n <et 2N far t 6 [0, i ]. 

Elementary arguments prove the existence of a constant ^4 > such that 

e^Ioit) <1- At forie[0,i ]. 

Therefore, 

to / oo \ *° 

I e t (e- t I (t)) d - 1 a » <2 " f ^ ecte '° / * 2iV " 1 ( 1 - ^) d_1 ^ < ee t0 ^/A, 

Vn=W / f )J 



so then 



0< Je-^Ioitf-Vj-Je-i* [(j2a n tA - 1 ] j < eeH^/A. (72) 



Observe that the upper bound in ([72"]) is independent of d. When combining (|7"Tj) with (f?2")) . we 
conclude that one can estimate 



|e- dt (/o(i) d -l)f by /e^n^a^-J - 1 J ^ (73) 
with error which is independent of d. 

It remains to study the asymptotic behavior in d of the second integral in fTB")) . For this, one 
needs to expand the <i-fold product of the polynomial integrand and carry out the integration. 
The integral of each term can then be expressed in terms of the incomplete Gamma function, 
which in turn can be estimated by the Gamma function itself. One then would employ elementary 
but somewhat involved expressions for multinomial coefficients, ultimately obtaining the following 
result. For any integer k > 1, there is a polynomial Pk(x) of degree k with Pfe(0) = such that 

j 'e- dt \(f^a n t 2n ^ -l ] jj=P k (l/d) + 0(d~ k - 1 ) asd^oc, 

We leave the details of these computations to the interested reader. 

7.3. The spectral zeta function. In very general circumstances, the spectral zeta function 
is defined as the Mellin transform of the thcta function formed with the non-zero eigenvalues. 
Specifically, for w £ C with Rc(ui) > d/2, we have that 

oo 

OvH = M9 N (w) = -L / (M*) - 1) t W T- ( 74 ) 
T(w) J t 
o 

As in section 5, let us consider a degenerating sequence of discrete tori. The elementary change 
of variables 1 1— ► u 2 t in (|74|) yields the expression 

oo 

aw^/ (e N {uh)-i)r^. (75) 
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Proposition 15.21 establishes the pointwise convergence of the integrand in ([75)) , and Lemma 
proves a uniform upper bound. Combining these results, as in the proof of Proposition 15. 4i 
together with ([75]). we conclude that 

lim (u- 2w ( N{u) (w)) =Ca(«0 (76) 
for any w £ C with Re(iy) > d/2 and where 

oo 

CA{W) = T^W) I (QaW-I)*' * (77) 
o 

is the spectral zeta function on the real torus AZ d \Z d . The difficulty is determining the cor- 
rect generalization of (|76p for w G C with Re(tu) > 0. The spectral zeta function ([75]) can be 
meromorphically continued to u> G C with Re(iu) > by writing 

oo 



T(w) 



i 



i 

f (0 N(u) {u 2 t) - V(N)e- 2d " 2t (I (2uH)) d ) t w j (78) 



T(w) 



l 



V(^)e^ 2<ill2 *(/o(2u 2 i)) d - l) t w j. 



T(w) 
o 

Similarly, the mcromorphic continuation of (|77|) is obtained by the expression 



Ca(w) 



r( 



oo 

w) J t 
l 



l 

f (e A (t)-V(A)(A7rtr d ^t w j (79) 





1 



l —J (v(A)(4nt.y d/2 -l) t w j. 



n 



The convergence of the first integral in ([78)) to the first integral in (|T9[) follows from Proposition 
5.21 and Lemma 15.31 as in the proof of ([76]) . The convergence of the second integral in ([75)) to the 
second integral in ([75]) follows from Proposition [531 using The Lebesgue Dominated Convergence 
Theorem. Let us write the third integral in ([78)) as 



(v(N)e-^ 2 %(2u 2 t) d l) r| = ™ y ( e - 2t / (2t)) - ± (80) 

o 

and the third integral in (|7T[) as 



^X^r" 3 - 1) r% = J 2 \ A) -- -. (si) 



Combining all of these results, we have shown that 
/ 



lim u~ 2w 



I WfrT \=U W )- V}*\ I ^ V (82) 



V 



T{w) J v v uv "' t ^ y ' (47r) d / 2 (w - d/2)T(w)' 
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Using Proposition 14.71 one can easily show that l|82p is equivalent to ([75]) when Re(ui) > d/2, 
noting that ([82]) holds for all w e C with Re(w) > provided Re(w) ^ d/2. By continuity, (|82|) 
extends to all w £ C with Re(w) > provided one uses the interpretation 

lk ?,o \ Ca{w) - u /2 /^l , F f \ ) = CT w=d/2 ( A (w), (83) 
\ (iiT) a i z yw — d/2)l (w) J 

the constant term in the Laurent expansion at w = d/2. 

The above computations in the special case d — 2 with w = 1 and w = 2 were obtained in section 
3 of [9]. In the case when w = 1 and d — 2, so then to = d/2, the results in [9] were expressed in 
terms of the modular forms. However, in light of the general Kroneckcr limit formula for Epstein 
zeta functions and its functional equation, as stated in section 12. 9( it is evident that for general 
d =/= 1, the limiting value obtained in (|83j) can be expressed in terms of a modular form (|22| . 



7.4. The Epstein-Hurwitz zeta function. Further analysis in [5] involves the investigation 
of the finite product Yi( s2 + Aj) f° r degenerating families of two-dimensional discrete tori for 
general s G C. If s — 0, then Theorem 15.81 determines the asymptotic behavior of the product of 
non-zero eigenvalues. For s / 0, we one includes the zero eigenvalues, which trivially introduces 
the multiplicative factor of s 2 . In this setting, we recall that Theorem 13.61 establishes a relation 
for J}(s 2 + Aj) in terms of integrals involving i-Bessel functions. Using the substitution t u 2 t, 
we arrive at the expression 

log~[[(( S /u) 2 + Aj) = ]T log((.s/ U ) 2 + Aj) - V(N/u)l d (s/u)+n N (s/u). (84) 

A j7 t0 



where 



and 



Us/u) = I (e- 2du2t e- s2t (I (2uh)) d - e^) j 



n N (s/u) = -l e~ sH (e N {uH) - V{N)e- 2duH {h{2uH)) d ) j - log( S 2 ). 
o 

Trivially, one has that 

Y[((s/u) 2 +A J )=u-^l[( s 2 + u 2 -A J ). 

We can now employ bounds from section 5 in order to determine the asymptotic behavior of I 
as u — > oo. Specifically, one uses Proposition 14. 7i Proposition 15. 21 Lemma 15.31 and Proposition 
15.51 to show that for any s £ C with Re(s 2 ) > we have the asymptotic formula 



~J{{s/u) 2 + Aj) = V{A/u)l d (s/u) + Ua{s) + o(l) as u -> oo 

where 

oo 

H A (s) = - I e- sH (64*) - V(A)(4nt)- d / 2 ) j - log( S 2 ). 
o 

The function H.a(s) is related to the regularized harmonic series, as defined and studied in |13| 
associated to the set {Aaj + s 2 } where {Aa,j} is the set of eigenvalues of the Laplacian on the 
real torus AZ d \R d . The general results from [13] establish that TLa(s) can be expressed as special 
values of the Epstein-Hurwitz zeta function formed with the set {Aaj + s 2 }. We refer the reader 
to [13] for further details and identities. 

7.5. General discrete tori. Let B be a positive definite d x d integer matrix and consider the 
discrete tori DTb = BZ d \Z d . The results from [15] easily extend to compute the spectrum of the 
Laplacian on DTb in terms of the dual lattice B* . The existence and uniqueness of the associated 
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heat kernel on DTb follows from general results (see, for example, [5] and [7]), thus allowing one to 
extend the results of section 3 above. Assume there is a one-parameter family DT B i u \ of discrete 
tori parameterized by u £ Z such that B(u) ju — > M as u — > oo where M is a positive definite d x (i 
matrix. The results of the present article apply when considering spectral invariants on DT B ^ u y 
For example, Theorem l5 . 81 will extend to obtain the asymptotic behavior of the determinant of the 
Laplacian on DTb( u ) with second-order term equal to the zeta-rcgularizcd spectral determinant on 
the real torus M Z d \R d . Additional results of the present paper, specifically the contents of section 
17.31 and section [7.41 carry through using the proofs given and only a slight change in notation. 

7.6. Height functions. The asymptotic expansion from Theorem 15.81 and its extension as 
outlined in section T7.51 can be interpreted as saying that logdet* is a height function on 
the space of discrete tori of fixed dimension since it tends to +00 through degeneration. In the 
case of real tori, one needs to introduce a minus sign and study — logdet* in order to have 

a height function which tends to +00 through degeneration; see, for example, [5]. The problem of 
finding real tori with minimum height remains a question of interest, see, for example, |23j . One 
point of future investigation is to see to what extent Theorem 15.81 allows for a connection between 
the problems considered in [23] and the study of the height function logdet* A.dt,* on discrete 
tori, which has the advantage of being an invariant which is defined as a finite product and not 
through meromorphic continuation. 
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